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CN ! ABSTRACT 
U 

^ l' Context. Because of new telescopes that will dramatically improve our knowledge of the interstellar medium, chemical models will 

' have to be used to simulate the chemistry of many regions with diverse properties. To make these models more robust, it is important 

^ , to understand their sensitivity to a variety of parameters. 

. Aims. In this article, we report a study of the sensitivity of a chemical model of a cold dense core, with homogeneous and time- 

independent physical conditions, to variations in the following parameters: initial chemical inventory, gas temperature and density, 
cosmic-ray ionization rate, chemical reaction rate coefficients, and elemental abundances. 
, Methods. We used a Monte Carlo method to randomly vary individual parameters and groups of parameters within realistic ranges. 

■^L ■ From the results of the parameter variations, we can quantify the sensitivity of the model to each parameter as a function of time. Our 

r H ' results can be used in principle with observations to constrain some parameters for different cold clouds. We also attempted to use the 

' Monte Carlo approach with all parameters varied collectively. 

f~| , Results. Within the parameter ranges studied, the most critical parameters turn out to be the reaction rate coefficients at times up to 

^ |i 4 X lO"" yr and elemental abundances at later times. At typical times of best agreement with observation, models are sensitive to both 

I ' of these parameters. The models are less sensitive to other parameters such as the gas density and temperature. 

O ' Conclusions. The improvement of models will require that the uncertainties in rate coefficients of important reactions be reduced. 

As the chemistry becomes better understood and more robust, it should be possible to use model sensitivities concerning other 



i-i 



, parameters, such as the elemental abundances and the cosmic ray ionization rate, to yield detailed information on cloud properties and 

^ , ' history. Nevertheless, at the current stage, we cannot determine the best values of all the parameters simultaneously based on purely 

observational constraints. 

Key words. Astrochemistry - Molecular processes - ISM: abundances - ISM: molecules - ISM: individual objects: L134N,TMC-1 
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1. Introduction sitivity of simple gas-phase chemical simulations of cold dense 
' interstellar cores to these parameters. The models are based on 
^ Molecules are powerful observational tools for the study of the ^j^^ pseudo-time-dependent approximation, in which the chem- 
O . physical conditions and dynamics of star-forming regions. Each ^^^^^ ^^^j^g^ ^^^^^ ^^^^ homogeneous density and tem- 
2 ■ step of the formation of a stellar and planetary system is char- ^^^^^^^^ Homogeneous models are also referred to as OD 
; acterized by a chemical composition, which directly reflects the (^ero-dimensional). Although more physically reasonable mod- 
> physical conditions of the medium and its evolutionary stage, ^j^^^j ^ j^^j^^^ formation of cold cores from more diffuse 
The far-infrared space telescope Herschel and (sub)-millimeter ^^^^^-^^ ^j^^ chemistry progresses, such models are rare given 
interferometer ALMA promise to open new windows on the ^j^^ j^^,^ unanimity of how the collapse proceeds. 
^ wealth of information found in the interstellar medium (ISM). ^^^^^ ^^e parameters in pseudo-time-dependent/OD mod- 
■ ■ • With the improvement of observational instrument sensitivity determined only to a limited extent by observations or their 
and resolution, and the opening of new wavelength ranges, more ^ompai-isons with model results are ( 1 ) the gas kinetic tempera- 
molecules will be detected in the ISM and many molecules will ^^^^ ^^^^i^y^ ^2) the cosmic-ray ionization rate, (3) the ele- 
be detected under a wider variety of physical conditions. The ^^^^^^ gas-phase abundances, and (4) the initial chemical inven- 
more molecules detected and the more diverse physical condi- ^^^^ ^j^^^^ parameters the rate coefficients (5) of the 
tions discovered, the more complex chemical models will have ^^^^ chemical reactions despite the fact that they are best de- 
to be to reproduce observations. termined by laboratory experiments rather than by observational 

The results of chemical models, or simulations, depend on constraints. We discuss the parameters in turn, 
a number of parameters or groups of parameters often poorly 

constrained. In this paper, we will be concerned with the sen- 1- The gas kinetic temperature and density are usually deter- 

mined by the excitation conditions of observed molecular 

* The osu.03.2008 network is available in electronic form at the emission lines. In addition to the uncertainties in the radia- 

CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via tive transfer analysis of the lines, inhomogeneities of the gas 

|http://cdsweb.u-strasbg.fr/cgi-bin/qc at?J/A+A/| may exist along the same line of sight. 
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2. There is no direct way to measure the cosmic-ray ionization 
rate ^ in dense clouds. There are two different approaches 
to constrain The first one, usually used, is to determine 
it by comparison of modeled and observed abundances for 
specific ionic species or their ratios. Currently, there ap- 
pears to be a dichotomy between diffuse clouds, in which 
the high abundance of H3 leads to a large value of and 
dense clouds, in which the ionization condit ions often lead 
to a value of 4' 1-2 orders of magnitude lower dWootten et al 



1982VCaseni et al.'.'l998^.'M cCall et al.Ll2003HLe Petit et al 
2004; Indriolo et al., 2009). There is some indication that 
the difference is simply due to the inability of low energy 
cosmic rays to penet rate into the center of dense clouds 
dPadovani et al.Ll2009l) . A second approach is to use the value 
measured in the solar system (I Spitzer & Tomaskol Il968t 
[Webber, 1998) and consider that f is constant in the Galaxy. 
However, in the solar system, the low energy end of the cos- 
mic ray spectrum cannot be determined directly because of 
the solar wind. 

The computed abundances of atomic and molecular species 
depend on the choice of elemental abundances in the gas 
phase. El ements are me asured in absorption in the diffuse 
medium ( Jenkinsll2004l) . In order to reproduce the observed 
gas-phase composition of dense clouds, it is often assumed 
that heavy atoms, including S, Si, Fe, Na and Mg, are de- 
pleted further from the gas between the diffuse (or translu- 
cent) and dense phases of a cloud life. Specifically, the ini- 
tial abundances of these elements in pseudo-time-dependent 
models are taken to be 2 to 3 orders of magnitude lower than 
in diffuse clouds but the efficiency of the depletions onto the 
grains is quite uncertain (Grae del et al., 1982). Moreover, the 
depletions of the elements may depend on the depth into the 
molecular cloud. 

Since the f = time of pseudo-time-dependent models is ar- 
tificial, the initial chemical inventory chosen can be artificial 
as well. Chemical models usually start with all the elements 
in the atomic form, except for H2, assuming that the den- 
sity of the gas had suddenly increased from a diffuse to a 
dense medium. Shock models of the conversion of diffuse to 
dense gas show, on the other hand, that the process is slow 
and that a significant amount of CO may be sy nthesized be- 
fore a sizable dense cloud core is produced (iBergin et al.L 
I2OO4I) . Preliminary results show that the formation of com- 
plex molecules is far slower if the starting form of carbon is 
CO (Hassel, Herbst, & Bergin, in preparation). Thus, the ini- 
tial chemical abundances and the cloud age, as determined 
by comparison of observational and model abundances, are 
correlated. 

Finally, the chemistry in the dense ISM involves thousands 
of gas-phase reactions, especially when one is dealing with 
complex molecules. Only a small fraction of these reactions 
have been studied at the low temperatures present in cold 
dense cores. Thus, even for systems studied in the laboratory 
at higher temperatures, there can be a difficulty in extrapolat- 
ing results to significantly lower temperatures. The most im- 
portant reactions at assorted cloud ages have been discussed 
in previo us papers with our sensitivity approach and those 
of others (Vasvunin et an.l2004tlmkelam et al.Ll2005ll2006t 
IVasvunin et al.. .20081)! 



In order to quantify the reliability of the chemical models 
and to determine how best to improve them, we need to know 
how sensitive they are to these parameters so that we can reduce 
the uncertainties in the more important ones. For example, if it is 



found that model results are more sensitive to current uncertain- 
ties in chemical rates than to elemental abundances, then future 
measurements of such rates should rank at high priority. If the 
opposite is true, priority might be given to developing better dy- 
namic models of clouds to follow the depletion of elements from 
the gas more appropriately. Moreover, if the models are com- 
pared with well-known cold clouds such as TMC-1 and L134N, 
the observations can be used to further constrain some of the 
model parameters within their computed uncertainties. In this 
article, we present sensitivity analyses of gas-phase cold dense 
cloud chemistry to all of the parameters discussed above. Details 
on the model are given in Section 2, while the sensitivity method 
and the ranges of values over which the parameters are varied 
are described in Section 3. The results of these analyses are dis- 
cussed in Section 4 while the constraints that observations of 
specific cold cores can yield concerning the model parameters 
are presented in Section 5. Finally, Section 6 contains a discus- 
sion. In our sensitivity calculations, one can proceed via the in- 
dependent variation of parameters or groups of parameters, or 
one can attempt to vary all parameters collectively. Although the 
latter approach is attempted, as discussed below, most of our use- 
ful results stem from the former treatment. 

2. Chemical model 

We used a new versio n of the Nahoon gas-phase chemical code 
(IWakelam et al.Ll2005h . which is significantly faster than the pre- 
vious version. The code is written in Fortran 90 and the differen- 
tial equation solver is DLSODES from the ODEPACK package 
(http://www.netlib.org/odepack/opkd-sum). Although the code 
is written for ID models (models that can be heterogeneous in 
one dimension), we used it in OD with gas temperature and to- 
tal hydrogen density set at the standard values for cold dense 
cores of T =10 K and «h= 2 x 10'* cm"-'. The visual extinc- 
tion is set to 10 so that photo-chemistry with external UV pho- 
tons is of little importance. The cosmic-ray ionization rate 
is 1.3 X 10"'^ s~'. We used the osu.03.2008 chemical network 
(http://www.physics.ohio-state.edu/~eric/research.html), which 
contains 13 elements, 455 species, and 4508 gas-phase reac- 
tions. This particular version of the osu network does not con- 
tain molecular anions, nor any dynamic depletion of gas-phase 
species. No reactions on grains are considered except for the pro- 
duction of H2 from two H atoms, a process that is treated in 
a pseudo-gas-phase manner by relating the density of grains to 
the overall gas density. The osu.03.2008 chemical network used 
in this work is available at the CDS with the following format. 
The first seven columns contain the reactants (three columns) 
and the products (four columns) of the reactions. The next three 
columns {a, [3 and y) give the parameters to compute the rate 
coefficients. Column 1 1 indicates the type of reaction. Column 
12 is the number of the reaction within a given type of reaction 
whereas column 1 3 is the general number of the reaction in the 
network. Finally, the last column gives the uncertainty factor for 
the rate coefficient. More details about the format can be found 
at http://www.physics.ohio-state.edu/~eric. Species are initially 
all in the atomic form except for H, which is 100% in H2. Table[T] 
summarizes our standard model. 

Listed in Table |2] the base elemental abundances used are 
oxygen -rich, and repres ent a low-sulfur version of the set labeled 
EA3 bv lWakelam & He rbst (2008). In this set of abundances, the 
amount of helium is based on the observed value in the Orion 
Nebula (Osterbrock et al., 1992). The oxygen, nitrogen and car- 
bon abundances seem to be rather constant in diffuse clouds 
and we have taken the values observed in ^ Oph dCardelh et al.L 
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Table 1. Standard dense cloud model. 



Parameter 


Value 


T 


10 K 




2 X 10'' cm-^ 


Ay 


10 




1.3 X 10-" S-' 


Initial abundances 


atomic except for 100% H2 


Rate coefficients 


osu.03.2008 



Table 2. Elemental abundances with respect to total hydrogen 
nuclei. 



Species 


Abundance 


He 


0.09 


N 


7.6 X 10"^ 





2.56 X 10-" 


F 


6.68 X 10-" 


C+ 


1.2 X 10-* 


S+ 


8.0 X 10-* 


Si+ 


0.0 


Fe+ 


1.5 X 10-* 


Na+ 


2.0 X 10-^ 


Mg+ 


0.0 


C1+ 


1.8 X 10-* 


P+ 


1.17 X 10-^ 



Table 3. Variational ranges of the parameters. 



Parameter 


Range 


N runs 


Reaction rate coefficients 


Uncertainty Factor 


2500 


Temperature* 


5-15 K 


2500 


H density* 


(1 - 3) X 10-* cm-3 




Elemental abundances 


+ 50% 


2000 


Cosmic-ray ionization rate 


(0.5 - 5.0) X 10-" S-' 


1000 


Initial Concentrations 


see text 


3500 



Varied together 



ll993HMever et al.Lll998h . Since we do not include gas depletion 
in Nahoon, we then removed 2.4% of C and 13% of O to ac- 
count for the fraction depleted on grain mantles in the f orm o f 
CO and H2O, as suggested by IShalabiea & GreenbergI (Il995h . 
For Na^, Cl^ and P ^, we use the observed values in ^ Oph 
d Savage & SembacR [1996). In dense clouds, iron, silicon and 
magnesium show additional depletion compared to the diffuse 
clouds, and we adopted the deple ted values recommended by 

Sulphur i s the m ost prob- 
Ruffle eiali Il999h . The low 
low- 



iFlower & Pineau des ForetsI (12003 



lematic element (see for instance 

value, taken from Graedel et al. (1982), is the standard 
metal" abundance. In summary, our low-sulfur EA3 abundances 
contain probable elemental abundances based on observational 
constraints. One of the consequences is that we do not produce 
large abundances of complex molecules since the C/0 ratio is 
significantly less than unity. These base abundances will be var- 
ied to determine the sensitivity of the model to elemental abun- 
dances. 



3. Sensitivity analyses 

3.1. Variations for most parameters 

The major goal of this analysis is to quantify the sensitivity of 
chemical models to variations of the parameters within realis- 
tic ranges of values. The method consists in generating new 
sets of parameters (P„) using random functions and comput- 
ing the corresponding chemical fractional abundances vs time 
(Xnit)). Using such an approach, we studied the sensitivities to 
the reaction rate coefficients, the gas temperature and density, 
the elemental abundances, and the cosmic -ray ionization rate. 
The variational range for each parameter is given in Table [3] 
For the reaction rate coefficients, we have assumed an ampli- 
tude variation equal to their uncertainties, as originally included 
in the UMIST database (see http://www.udfa.net/) and subse- 
quently in the osu database , with a log-normal distribution (see 
IWakelam et all. 120051 12OO6I) . We define Fi, the uncertainty factor 
of the rate coefficient, at a Icr level of confidence. As a conse- 
quence, there is a 68.3% chance that the rate coefficient lies be- 
tween ki/Fi and fe; x Fi, where ki is the listed rate coefficient for 
reaction /. For the other parameters, except for the initial chemi- 
cal inventory (see below), flat distributions (no preferred values) 
were assumed with ranges determined by reasonable uncertain- 
ties, which are mainly +50% from the center or peak value. The 
gas temperature then lies between 5 and 15 K, the total hydro- 
gen density between 1x10'* and 3 x lO'* cm--', and the elemental 
abundances between +50% of their standard values. The cosmic 
ray ionization rate, which is clearly uncertain to more than a fac- 
tor of two, is varied between 5 x 10-'** and 5 x 10-'^ s-' . 

The variations of the elemental abundances are sufficient to 
probe perturbations around the chosen values, including C/O ra- 
tios greater than unity, but not large enough to include the to- 
tal diversity used by previous modelers. All individual parame- 
ters (e.g., the cosmic ray ionization rate) or groups of parame- 
ters (e.g., reaction rate coefficients, elemental abundances) were 
varied one at a time except for the temperature and the density, 
which were varied together although they are not correlated. In 
total, we ran 8 000 different models of the dense cloud in vary- 
ing these parameters. We used a Monte Carlo method to sample 
parameters randomly; this is most efficient for groups of param- 
eters such as the reaction rate coefficients. For individual param- 
eters such as the cosmic-ray ionization rate, a non-random sam- 
pling method could have been used to span the range of values. 

To quantify the dispersion of the computed abundance pro- 
file of a specific species at a given time when one parameter 
or a number of different parameters is varied, we calculate the 
standard deviation of the profile in logarithmic units using the 
formula 



(r(f) = 



(1) 



where the sum is over the number of computed abundances A^, 
X„(t) is the abundance of the specific species for run « at time 
t and logiXit)) is the mean of the logarithm of the abundance 
at time t. We have chosen to compute cr in logarithmic units in 
order to deal with numbers over many orders of magnitude. In 
addition, the uncertainties of most rate coefficients are given as 
factors rather than additive terms, making logarithmic units more 
natural. If the abundance profile at time t follows a log-normal 
distribution (Gaussian in logarithmic units) with cr(f) - 1,68.3% 
of the abundances lie within a factor of 10 from the mean abun- 
dance. In the rest of the paper, we will use this definition of the 
standard deviation. 
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Fig. 1. CO abundance computed in a diffuse/translucent cloud. 
Left panel: plot of 10% of the runs. Middle panel: the gray lev- 
els represents the percentage of runs. Any level of gray: at least 
one run. Two darkest gray levels: 95% of all runs at one given 
time (2 sigma if the distribution is log-normal). Darkest gray: 
68% of all runs (1 sigma if the distribution is log-normal). Right 
panel: distribution of the fractional abundance at 10^ yr, where 
the median value is approximately 2 x 10"^. 



3.2. Determination of the initial ctiemical inventory 

The initial atomic and molecular abundances for the cold core 
are treated uniquely. To study the sensitivity to the initial concen- 
trations, we need to have realistic distributions for dense cloud 
models. Since dense molecular clouds are presumably formed 
from diffuse or translucent gas, we first compute the steady-state 
chemical composition of such a gas with a temperature of 50 K, 
a total hydrogen density of 2 x 10^ cm"^ and a visual extinc- 
tion of unity. The species start in their atomic form prior to the 
diffuse/translucent cloud stage, except for hydrogen, which is as- 
sumed to be divided equally among atomic and molecular forms. 
The fixed elemental abundances used for the diffuse/translucent 
stage are the ones listed in Table |2l although a more realistic set 
would have less depletion. In order to get realistic distributions, 
we randomly varied the reaction rate coefficients within their un- 
certainty ranges using an additional 3500 runs. The abundances 
and distributions we obtained at steady-state were then used as 
the initial chemical inventory for the dense cloud. As an ex- 
ample, the CO abundance distribution in the diffuse/translucent 
cloud conditions is shown in Fig.[T](left and middle panels) as a 
function of time. The distribution of CO at steady-state (10^ yr) 
used as the initial distribution for the dense cloud is shown on 
the right panel and can be fitted by a log-normal distribution. It 
peaks at a fractional abundance of ^ 2 x 10"^ with a standard 
deviation of a factor of 6.3. For comparison, the C^ and C dis- 
tributions have their steady-state peaks at 6 x 10"^ and 10"^ 
respectively, so that the gas is still mainly atomic except for H2, 
which takes up most of the elemental hydrogen; the atomic H 
distribution peaks at a; 3 x 10~^. These results are somewhat dif- 
ferent from our typical assumption that carbon and oxygen are 
purely atomic and hydrogen purely molecular at the initial stage 
of a dense cloud. Our distributions span two extreme scenarios: 
one in which most of the carbon is still in the atomic form (C^ 
abundance 1 x 10 '*) and one in which CO is the dominant 
form of carbon (CO abundance ^ I x 10 '*). These extreme val- 
ues can arise depending up on parameters in the shock models 
for dense cloud formation of iBergin et al.1 (|2004|) . 



3.3. Variation of all parameters 

In the sections presented above, we have emphasized how we 
vary the parameters or groups of parameters one after the other 
to study their eff'ect independently. The alternative approach is to 
develop a model in which all parameters are varied collectively. 
To achieve this goal, we have once again used a two-step pro- 
cedure involving the diff'use and dense cloud stages, discussed 



Step 1 : Diffuse cloud 

Fix parameters: T, n(H) and C, 
N sets of elemental abundances 
N sets of rate coefficients 

N sets of steady-state abundances / 



Step 2: Dense cloud 

N sets of T, n(H) and C, 
N sets of rate coefficients 
N sets of initial chemical 
/ composition 

\ 

N chemical evolutions 



Fig. 2. Steps used to vaiy all the parameters of the model collec- 
tively. 



below. The procedure is also summarized in Fig.|2]for better un- 
derstanding. 

For the first step, we once again compute the steady-state 
chemical composition of a diffuse cloud with a fixed tempera- 
ture of 50 K, a density of 2 x lO"* cm"-', visual extinction of 
unity, and cosmic-ray ionization rate of 1.3 x 10"'^ s"' (as de- 
scribed in the previous subsection). But now we vary at the same 
time the elemental abundances and the rate coefficients within 
the ranges given in section [TTl using the Monte-Carlo approach 
with a different distribution of random numbers for both sets of 
parameters. We obtain - 3500 different steady-state chemical 
compositions as a consequence of varying the sets of elemental 
abundances and rate coefficients simultaneously. In the second 
step, we use this distribution of 3500 chemical compositions as 
the initial inventory and compute the chemical evolution of a 
dense cloud using the previous 3500 sets of rate coefficients. In 
addition, we vary simultaneously the temperature, the density 
and the cosmic -ray ionization rate using different random num- 
ber distributions for each parameter We then obtain 3500 dif- 
ferent chemical evolutionary paths, each of which corresponds 
to to a different temperature, density, cosmic-ray ionization rate, 
set of rate coefficients, set of elemental abundances and initial 
chemical composition. 



4. Results 

In subsections 4.1-4.4 below, we discuss the sensitivities of the 
abundances to individual variations of parameters or groups of 
parameters with the others set equal to their standard values. 
Afterwards we compare these sensitivities and also consider the 
sensitivities to all parameters varied collectively. 

4.1. Sensitivity to the initial concentrations and their 
distributions 

The sensitivities of selected dense cloud abundances to changes 
in initial chemical concentrations are shown in Fig. [3] as func- 
tions of time. The initial distributions of the molecular concen- 
trations introduce spreads in the time-dependent solutions that 
last until the system is at steady state, where abundances no 
longer depend on the initial concentrations because we do not 
have any bistability for the parameters used. Although some dis- 
tributions, such as that for atomic H, show standard deviations 
(cr = 0.23, see section ITTI ) as large as a factor of ^ 2 for times in 
the dense cloud stage up to 10^ yr, most distributions are more 
sharply peaked, as can be seen in Fig. |3] For the species shown 
in this figure, the actual logarithmic standard deviations cr for 
O2, CH, C3O and HC7N are 0.14, 0.12, 0.07 and 0.17 respec- 
tively at 10^ yr For comparison, we have started the dense cloud 
stage with our standard initial chemical abundances - atoms ex- 
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2345678 2345678 



log(time[yr]) log(time[yr]) 

Fig. 3. O2, CH, C3O and HC7N abundances as a function of time 
in the dense cloud stage, starting with a distribution of atomic 
and molecular abundances discussed in the text. The black lines 
represent abundances obtained if we use our standard atomic ini- 
tial concentrations except for H2. 




0.01 0.1 1 10 100 1000 -5.5 -5.0 -4.5 -4.0 -3.5 -3.0 



Fig. 4. Left panel: Calculated HC3N, C5H and CH3OH abun- 
dances as a function of the initial concentration ratio C^/CO of 
the dense cloud stage. Right panel: Calculated CH2 and HC3N 
abundances as a function of initial atomic hydrogen fractional 
abundance. Results are shown for 10^ yr. 



cept for 100% H2 - at the elemental abundances found in Table|2l 
The fractional abundances obtained with these initial concentra- 
tions are also plotted in Fig. [3] with single black lines. For most 
molecules, there is very little difference between the black lines 
and the distributions at any time, although what differences exist 
occur mainly at earlier times. For complex molecules however, 
such as HC7N, CH3OH, and NH2CHO, which typically peak at 
so-called early times in the range lO"*"^ yr, the early-time abun- 
dances are significantly higher with the standard than with the 
calculated initial abundances. 

One reason may lie in the relatively high abundance of 
atomic hydrogen calculated at steady-state in the diffuse stage. 
Unlike the case for standard abundances, in which, starting from 
zero, the fractional abundance of atomic H reaches a constant 
10"* by 10^ yr, the calculated atomic hydrogen distribution at 
steady-state in the diffuse stage remains peaked at * 3 x lO""' for 
up to 10^ yr in the dense cloud stage. Note that our model in- 
cludes the H2 an d CO self-shielding using approximations from 
iLee et al ] dl996h . The excess atomic hydrogen calculated in the 
diffuse cloud stage and used initially in the dense cloud stage can 
interfere with the production of complex molecules by depleting 
some intermediate hydrocarbon radicals via neutral-neutral reac- 
tions. Another possible reason is the relatively high abundance 
of CO calculated to occur at steady-state in the diffuse stage. 



Once a significant amount of carbon becomes entrained in CO, 
it is not as easy to form complex carbon-containing species. 

To determine which of these effects is the more important 
for the dense cloud stage, we first varied the initial C^/CO dense 
cloud ratio over a wide range while fixing the overall carbon 
abundance as well as the other elemental abundances at the stan- 
dard values shown in Table |2] We separately did an analogous 
variation for the initial H/H2 ratio in the dense cloud stage, in 
which case all the carbon lies initially in C^. Some results at a 
dense cloud time of 10^ yr can be seen in Fig.|4] Larger carbon- 
containing organic species are clearly enhanced for larger initial 
C^/CO ratios, but there is a saturation point for C^/CO larger 
than 1 . Because our standard initial C^/CO ratio exceeds unity 
and our calculated dense cloud initial abundance ratio has a me- 
dian value near 5, the variation of this ratio between our two 
sets of initial concentrations is not the reason for differences 
in abundances of complex molecules such as HC7N. Regarding 
the initial fractional abundance of atomic hydrogen, however, 
we see that even the smaller organic species CH2 and HC3N 
are reduced in abundance as the ratio increases to 10"^; the re- 
duction in abundance shown is even larger for more complex 
species. Note that there is also a saturation effect here. In particu- 
lar, the molecular abundances are not sensitive to H initial abun- 
dances smaller than 10"^, but this value is well below the range 
of values in our calculations. Within the ranges of initial con- 
centrations for the dense cloud stage studied, the major source 
of sensitivity is thus the residual atomic hydrogen abundance. 
For initial concentrations, however, in which there is more CO 
than C^, the abundances of larger organic species will be signif- 
icantly reduced at early times. The fact that a large abundance of 
atomic hydrogen can significantly reduce calculated peak abun- 
dances of complex molecules at early time is of some impor- 
tance considering that 21 cm absorption studies indicate signifi- 
cant abundances of atomic H in some c old cores that are higher 
than calculated by our standard model dRrco et al 112009). This 
possible discrepancy argues against the high calculated early- 
time abundances for complex molecules in the simple pseudo- 
time-dependent model and for its possible replacement by mod- 
els that take into account cloud formation and heterogeneity. On 
the other hand, the presence of anions or PAH's, not included 
in this study, tends to aid the synthesis of co mplex molecules 
(IWakelam & Herbsil2008HWalsh et alil2009l) . 

In the variations of the other parameters discussed below, 
we return to our standard initial chemical abundances for dense 
clouds. This choice maximizes complex molecule abundances at 
early times. Nevertheless, the abundances achieved are still low 
compared with low-metal e lemental abundance s, as can be seen 
in the recent calculations of IWalsh et al.l (l2009l) . 

4.2. Sensitivity to fhe reaction rate coefficients, temperature, 
and density 

The sensitivities of atomic and molecular abundances in cold 
dense clouds to variations in reaction rate coefficients and 
gas temperature and density have already been studied by 
Wakelam et alJ (I2006h . For this reason, we will just mention 
the salient results of our new calculations. For rate coefficients, 
the computed abundances have log-normal distributions result- 
ing from the probab ility distribution func tions of the reaction 
rate coefficients (see IWakelam et aO, l2006l) . Moreover, a varia- 
tion of temperature within the range 5 - 15 K mainly affects the 
nitrogen-bearing species. In fact, much of the nitrogen chemistry 
starts with the reaction N^ +H2 — > NH^ + H, which is either 
slightly endothermic or possesses a small barrier. The rate coef- 
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Fig. 6. HC3N sensitivity to the cosmic-ray ionization rate (^). 
Left panel: Probability density as a function of time for the dif- 
ferent values of ^. Right panel: HC3N abundance as a function 
of ^ for three different times. 



ficient of this reaction at temperatures under 100 K depends upon 
the fraction of molecular hydrogen i n its excited ortho state; the 
results of more detailed calculations ( Le Bourlo ?. 1991) yield an 
effective activation energy of 85 K for the temperature range 
10-100 K, as listed in the osu network. Temperatures below 
10 K then produce smaller abundances of N-be aring spec i es tha n 
temperatures above 10 K (see Fig. 4 of Wake lam et all l2006h . 
Within the chosen range, the density has little influence on the 
abundances. 



4.3. Sensitivity to tlie cosmic-ray ionization rate 

Most of the molecules are sensitive to the cosmic-ray ioniza- 
tion rate (() since cosmic ray ionization is the main source of 
atomic and molecular ions for cold clouds, and the abundances 
of neutral species are determined mainly by reactions involving 
ions and/or electrons. Sensitive species can be divided into two 
groups: 1) those sensitive to a change in ^ at all times and 2) 
those only sensitive to f between 10'* and 10^ yr. The first group 
is composed of molecules that are more or less direct products 
of the ionization of H2 by cosmic rays, such as HJ, OH (see 
Fig.jSj and O2. Complex molecules, such as the cyanopolyynes, 
are part of the second group, as can be seen in Fig.|6]for the case 
of HC3N. Some molecular abundances increase as ( increases 
(for instance OH, CO, O2, C3, C3O, SO, NO, and most of the 
ions), whereas others clearly decrease with increasing ^ (for in- 
stance the cyanopolyynes, H2O, OCS, CS, CH3OH and HCN). 
The cases of OH and HC3N can be seen in the right panels of 
Figs. |5] and |6] Presumably, the distinction for neutral species 
arises from whether or not ions are more important in the for- 
mation or the destruction of molecular species. The species OH 
and H2O for instance are both produced by the dissociative re- 



-10 - 



-12 




7 8 



bg (time [year]) 



200 










A 


1 100 
E 






z 

50 











r . . 1 . 


20 


-15 -10 -5 
log[X(i)] 


-8 
q7 -10 


■ 11 




abundar 






logi 








0.1 


1.0 

Elemental C/0 



Fig. 7. HC3N sensitivity to elemental abundances. Left panel: 
HC3N probability density as a function of time. Right upper 
panel: distribution of HC3N abundance at 10^ yr. Right lower 
panel: HC3N abundance as a function of elemental C/O ratio at 
10^ yr 



combination of H^O^ but H2O is mainly destroyed by an ion- 
molecular reaction with C^ whereas OH is destroyed by reaction 
with neutral N. 

The de pendence on 4" at steady s tate has been d i scusse d pre- 
viously by Lepp & Dalgarno' (1996') and ' Lee et all (Il998l) . who 
found that when photodestruction by external photons is ignor- 
able, then fractional abundances depend on the ratio ^/hh- Thus 
the fact that the distribution for HC3N collapses to near a point 
at steady state despite the large variation in f shows that there 
is also very little dependence on riu whereas the distribution for 
OH shows that this radical has some (inverse) dependence on rin. 



4.4. Sensitivity to ttie eiementai abundances 

The effects of variations in the elemental abundances are shown 
in Fig. Q for HC3N, where the distribution spreads out with in- 
creasing time, reaching a standard deviation that exceeds 1.0, 
which implies a 68% chance of finding the abundance to be 
spread over 2.5 orders of magnitude, when steady-state is finally 
reached. The most important variations in the range considered 
are those for the elements carbon and oxygen, and most of the 
sensitivity is to the ratio of these elements, which varies between 
0.2 and 1.4. The correlation between the HC3N abundance at 
10^ yr and the C/O elemental ratio is shown in the right lower 
panel of the figure, where it can be seen that HC3N increases in 
abundance quite dramatically with C/O whereas the spread at a 
given C/O ratio caused by different values of C and O as well as 
variations of the other elemental abundances is smaller. In gen- 
eral, larger abundances of C-rich molecules such as HC3N are 
obtained with increasing C/O. The reason is that the greater the 
C/O ratio, the greater the amount of carbon that is not taken up 
in the form of CO. When the C/O elemental ratio becomes larger 
than 1 , we go from a chemical regime dominated by oxygen to 
a regime where carbon chemistry is the more important. Since it 
takes some time for a large amount of CO to be synthesized, the 
largest effect is seen at the longest times. 



4.5. Comparison of the different sensitivities 

To compare the effect of each parameter, including the initial 
concentrations as described in Section l4n we computed the log- 
arithmic standard deviation cr, of the abundance distributions 
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Fig. 9. Abundance of HC3N as a function of time for the case 
of all parameters varied together on the left and for the variation 
of rate coefficient only on the right. In each panel, the middle 
curve represents the median abundance while the two external 
curves represent a Icr logarithmic standard deviation from this 
abundance. 



at each time for each parameter group / - initial concentra- 
tions, elemental abundances, cosmic -ray ionization rate, temper- 
ature/density, and reaction rate coefficients - as well as for the 
case where we vary all parameters together. Note again that a 
logarithmic standard deviation of 1.0 refers to a factor of 10 in 
either direction (see Section lTTI ). As an illustration, we show the 
individual logarithmic standard deviations for several molecules 
as a function of time in Fig. [8] Different molecules are sensitive 
to different parameters and this sensitivity varies with time. For 
HC3N, the reaction rate coefficients dominate until late times 
when the species is more sensitive to the elemental abundances 
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(see Section |44] |. The species NH3, on the other hand, is more 
sensitive to the gas temperature and density for times larger than 
10^ yr (see Section 14.2b . As expected, the standard deviation 
when varying all parameters is generally larger than the ones 
for individual parameters but much smaller than the sum of the 
individual effects. 

In some cases, however, the standard deviation for all pa- 
rameters is smaller than for one of the individual parameters or 
parameter groups, such as the rate-coefficient ("R. C") case for 
HC3N between 10^ and 10^ yr This situation can happen be- 
cause the amplification of the error affecting the rate coefficients 
and its propagation into a spread of abundances depends on the 
set of parameters other than the rate coefficients. In some cases, 
the standard set of parameters will produce a dispersion that is 
larger than what is obtained with most of the other sets that are 
generated in the case for all parameters. As a consequence, the 
value of cr(f) obtained by varying only the rate coefficients can 
be higher than what is found when varying all the parameters 
together, which attenuates the dispersion of abundances. Large 
values of cr(f) are often associated with a stiff variation of the 
abundance, because one of the effects of changing the rate co- 
efficients is to produce a time shifting of the chemical evolu- 
tion. If a species experiences a stiff variation of its abundance at 
a time f, this time-shifting produced by varying the rate coeffi- 
cients results in a large dispersion of its abundances at t. When 
the standard set of parameters produces the stiffest variation, the 
deviation in the R. C. case can exceed that for all of the param- 
eters. In Fig. |9l we show the mean abundance of HC3N with 
its (J standard deviation in the case of the variation of all pa- 
rameters and in the case of the variation of the rate coefficients 
only. For each case, we show the standard deviation obtained at 
10"* yr by a gray vertical line. In the R.C case, the abundance of 
HC3N increases much more stiffly than in the case for all param- 
eters and this introduces a larger standard deviation at a specific 
time. In this specific case, the change in stiffness, which also 
produces a difference of the mean abundance, arises because of 
the sensitivity of the stiffness to the initial conditions. As dis- 
cussed in Section l4n the case of all parameters uses computed 
initial compositions whereas the R. C. case assumes the standard 
atomic composition, which results in a stiffer evolution around 
10"* yr In the case of NH3, as shown in Fig.|8] at 10^ yr the sen- 
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sitivity to the temperature produces a slightly larger cr than the 
case for all parameters at a time when the NH3 abundance in- 
creases stiffly from 10"^ to 4 x 10"^ on a very short period of 
time. At that time and for the particular set of parameters, the 
stiffness of the NH3 abundance curve is diluted by the variation 
of the other parameters. 

To study the general effect on the complete model, we cal- 
culated the average individual logarithmic standard deviations 
for all species as a function of time through 10^ yr, as can be 
seen in Fig. [10] Although the relative importance of the param- 
eters varies with time, sensitivity to the reaction rate coefficients 
dominates for ages less than 4x10^ yr For later times, the el- 
emental abundances are the main source of uncertainty in the 
model. At early time, where model results are typically closest 
to observed values for cold cores, the average logarithmic stan- 
dard deviations of the various parameters range from 0.5 (factor 
of three; reaction rate coefficients) to 0.1 (factor of 1.3; initial 
concentrations). Of course, these standard deviations are based 
on the ranges of the parameters. Although most of these ranges 
are reasonable, based on a variety of criteria, the range for the el- 
emental abundances does not cover some of the large differences 
between so-called high-metal and low-metal abundances. With 
larger ranges, the sensitivity to elemental abundances is likely to 
increase in importance even at shorter times. Another effect not 
included in the figure is that of the large difference between the 
range of initial H atom concentrations when the initial concen- 
trations are varied, and the standard initial abundance of zero. 
Even the standard deviation introduced by a variation of initial 
H between (standard initial abundance) and lO^^* is however 
smaller than all the standard deviations discussed in this section, 
because its main effect is on the most complex molecules. 

5. Comparison with observations 

It is common to use chemical models to constrain some of the 
physical parameters of a cloud (cosmic-ray ionization rate, age, 
etc.) by comparing observed abundances with models in which 
the selected parameter is changed. Note that the age of a cloud 
is not one of the parameters that is handled by our random treat- 
ment; it is possible but not necessary to handle time randomly. 
In any case, this method of constraint for one parameter or group 
of parameters requires that the other parameters be well known 
and/or that the model results, such as abundances, do not de- 
pend on the adopted values and uncertainties for the other pa- 
rameters. Both hypotheses are usually not completely true, es- 
pecially for reaction rate coefficients. Using our sensitivity anal- 
ysis, we have previously studied the possibility of putting some 
constraints on some of the model parameters for dense clouds by 
comparing our modeling results with observations in two dense 
clouds L134N, with 4 2 species, and TMC-1, with 53 species 
jWakelam et al.Ll2006l) . Our major approach was to define agree- 
ment by the number of species for which the calculated distri- 
bution of abundances overlaps with the observational value and 
uncertainty. 

To compare observed and modeled abundances here, we 
have instead computed at each time a function D{t), defined by 
the equation 

D(t)^Yj\^og{xljt))-\og{xij\ (2) 

with X^^^^ the observed fractional abundance of species j and 
X'^^^if) the abundance of species j computed by the model at 



time t. For each parameter or group of parameters, we then at- 
tempted to determine the minimum value of D{t) as a function 
of both the parameter and time. The smaller the value of D, the 
better the agreement. Our definition of D makes the "agreement" 
with minor species such as cyanopolyynes as important as the 
agreement with major species such as CO. There is, in our opin- 
ion, no reason why CO should have more weight in the compar- 
ison than minor species. Indeed, reproducing the abundance of 
a minor species might indicate that the network works correctly 
at a more detailed level. We could introduce a weighting factor 
if we had some indication that some observed abundances are 
wrong or highly uncertain for instance. This is generally not the 
case here. Nevertheless, it is indeed true, as the referee hints, that 
the parameter D is somewhat arbitrary. 

5. 1 . Constraints on the model parameters using all observed 
species 

Here we sum over all observed species to obtain D for a given 
source. An important result of our analysis is that the constraints 
on one specific parameter depend on the values of the other pa- 
rameters. This result is illustrated by Fig.[TT] which shows three 
separate plots. On the left one, the minimum value of D{t) for 
TMC-1 is displayed as a function of the elemental abundance 
ratio C/O for the case where we vary only the elemental abun- 
dances. The middle plot represents the same quantity but in the 
case where we vary all parameters together. As can be seen on 
this figure, there is a large spread for the minimum D in the 
middle panel, and to observe the trend better, we show on the 
right plot, the mean values of minimum D obtained in this case, 
as well as the minimum values as a function of elemental C/O. 
Although we see the same trend in both cases that the observa- 
tions in TMC- 1 are better reproduced by larger values of C/O, we 
also see that the best agreement (D « 45) also exists for smaller 
values of elemental C/O in the case where the other parameters 
are varied. In other words, the "best model" is not necessarily 
the one at the extreme value of a parameter. On the plots of the 
right of Fig.[TT]for instance, the best agreement is rather flat from 
C/O=0.4 to 1.4. We can see this effect in other instances. For ex- 
ample, if we only consider the agreement with the observations 
as a function of the initial ratio C^/CO while keeping the ele- 
mental abundances fixed and all the other elements in the atomic 
form, we find that observations in TMC-1 are best reproduced 
if all the carbon lies in the atomic form. If we change other ini- 
tial abundances, such as the initial H/H2 and/or the initial N/N2 
abundance ratios, we find the best agreement for other values of 
the initial C+/CO. 

We will not give here the set of parameters for the "best 
model" because we did not explore the full parameter space. 
Doing so would be too time-consuming. The best approach for 
the moment is to reduce the range of possible variation of the 
parameters, especially for the reaction rate coefficients, which 
seem to be dominant for many molecules. 

5.2. Constraints on the model parameters using selected 
species 

In the previous section, we discussed some limitations on the use 
of the abundances of all species observed in a cloud to constrain 
some of the model parameters via the method of minimization 
of D. One may do better by only using species sensitive to a 
given parameter so that their sensitivity is not masked by less 
sensitive species. Molecular ions, such as H^, have been used to 
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Fig. 11. Minimum D in TMC-1 as a function of the elemental ra- 
tio C/O. The left panel shows the results when we vary only the 
elemental abundances of the model, the middle panel shows the 
result when we vary all the parameters together, and the right 
panel shows the mean and minimum values from the middle 
panel. 



constrain th e cosmic-ray i onization rate in dif f use re gions (see 
for instance iMcCall et all 120031: iLe Petit et all l2004 . In dense 
clouds, the abundance ratios HCO"'/CO and DCO+/HCO+ have 
been proposed to constrain the cosmic ray ionization rate an d the 
fractional ionization (see for instance ICaselli et all Il998l) . We 
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Fig. 12. HCO^/CO abundance ratio as a function of ^ at 10^ yr 
in the case where we vary all parameters together. 



show in Fig.[T2|the abundance ratio HCO^/CO as a function of 
the cosmic-ray ionization ^ in the case where all parameters are 
varied at the same time. This ratio depends on time and Fig.[T2l 
is for 10^ yr, the so-called early time of best overall agreement 
by number of species. There is a clear correlation of the abun- 
dance ratio with ^ but the spread of possible HCO^/CO values 
at a spe cific i s at least o ne order of magnitu de. According 
to Dicke ns et al] (|2000) and lOhishi et all (|1992|) . the observed 
value of HCO+/CO is about 1.5 x lO^^ in L134N and lO""* in 
TMC-1 (CP). Compared with Fig. [121 these results mean that 
^ should be > 5 X 10 s ' in both clouds. The authors do 
not however give the uncertainties in the observed abundances. 
Uncertainties should indeed be given to improve the utility of 
sensitivity studies. The analysis of iCaselli et alj d 19981) also in- 
dicates a rather high value of ( in TMC-1 (CP) in the range 
(4 - 8) X 10 s 1. If we plot the HCO+/CO abundance ratio 
vs ^/n, a parameter discussed earlier, we obtain that (/n should 



exceed 5x10 
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6. Discussion 

We have studied the sensitivity of the chemistry of cold dense 
cores to an assortment of parameters using the OD gas-phase 
model Nahoon (Wakelam et al., 2005) in which chemical abun- 
dances evolve from given initial abundances. As variable pa- 
rameters, we considered the temperature and density of the gas, 
the cosmic-ray ionization rate, the elemental abundances, the 
reaction rate coefficients, and initial concentrations. These pa- 
rameters or groups of parameters were all varied within realis- 
tic ranges of values eith er individually o r simultaneousl y using 
Monte Carlo methods ( ^svunin et all 12004; Wakelam et all 
I2005L 120061: IVasvunin et all l2008h used previously mainly to 
study the sensitivity of molecular abundances to chemical reac- 
tion rates. For the case of the elemental abundances, our vari- 
ations were carried out randomly f rom a base of abundances 
based on a new set first defined in IWakelam & Herbsll (120081) 
as appropriate for cold dense clouds. Our results depend to some 
extent upon whether individual parameters are varied while oth- 
ers are held fixed at standard values or whether all parameters 
are varied simultaneously. Nevertheless, certain trends can be 
found. Among the most important is the finding that, when av- 
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eraged over all molecules, the dominant source of uncertainty in 
predicted molecular abundances for times less than 4 x 10^ yr 
is the uncertainty in rate coefficients, whereas the dominant un- 
certainty at later times is caused by uncertainties in elemental 
abundances. The former source of uncertainty can be reduced by 
improved laboratory and theoretical determinations of rate coef- 
ficients at appropriate temperatures, especially when sensitivity 
analyses point out the spec i fic rea ctions of greatest importance, 
as done in Wake lam et al. I (l2009h . The latter source of uncer- 
tainty is astrophysical in origin, and suggests that more attention 
be paid to a physical understanding of how gas-phase elemental 
abundances evolve as dense clouds are formed from more diffuse 
gas. 

Within the ranges shown in Table [3] the sensitivity to the 
other parameters is generally less. As long as the temperature 
and density of the gas uncertainties lie within their listed ranges, 
the molecular abundances predicted by the gas-phase model are 
quite robust; i.e., they do not depend much on T and hh- The 
only exceptions are the N-bearing species containing only nitro- 
gen and hydrogen, such as NH3, which are produced much less 
efficiently at temperatures lower than the standard value of 10 K. 
Many molecular abundances depend on the value of the cosmic- 
ray ionization rate, either at all times or only between cloud ages 
of 10"* and 10^ yr An example of the former is OH and one of the 
latter is HC3N, as shown in Fig.|6] Of all the parameters, the ini- 
tial concentrations with the chosen elemental abundances seem 
to be the parameter to which model results at all reasonable times 
are least sensitive. Nevertheless, we found that cloud chemistry 
has to start with a significant fraction of the carbon in the atomic 
form and a very small amount of the hydrogen in atomic form to 
synthesize complex molecules. Our standard model meets these 
requirements, but it may be in co nflict with with some recent 21 
cm studies by iKrco et aTl (l2009h . which show a larger amount 
of atomic hydrogen. The range of atomic hydrogen abundances 
used when the concentrations are varied differs strongly from 
the standard model in which all hydrogen starts in its molecular 
form. It is thus critical to understand both the physics and chem- 
istry of the stages of collapse leading to the formation of cold 
cores for this problem as well as the elemental abundance prob- 
lem. A new treatment of the gas-grain chemistry of cold cores in 
the process of formation is being prepared by Hassel, Herbst, & 
Bergin. 

In addition to the study of model sensitivity to assorted pa- 
rameters, we have utilized the simultaneous variations of these 
parameters to attempt to determine their optimum values by 
comparison of calculated and observational results for the cold 
cloud cores TMC-1 and L134N. In this paper, we used an indi- 
cator of agreement that minimizes the sum of the absolute val- 
ues of the difference of the logarithms of the calculated and ob- 
served abundances. Constraints on individual model parameters 
are highly sensitive to the values of the other parameters (espe- 
cially reaction rate coefficients) when variations are run simulta- 
neously. It is currently not possible to constrain all the parame- 
ters by comparison with observed abundances by rigorous meth- 
ods so that efforts have to be made to reduce the uncertainties 
for some of the parameters. Uncertainties in rate coefficients for 
instance can be reduced by laboratory measurements or theoret- 
ical calculations on reactions clearly identified as quantitatively 
important for the model predictions. 

Finally, although the strength of sensitivity methods for gas- 
phase chemical models has been clearly demonstrated in this and 
earlier studies for both cold clouds and more complex objects 
such as protoplanetary disks dVasvunin et all |2008), the meth- 
ods will need to be generalized to estimate the sensitivity to both 



chemical processes on grain surfaces and adsorption and desorp- 
tion processes. The use of sensitivity methods for surface chem- 
istry may become quite feasible within the next decade given 
the rapid pace of advance of both theory and experiments in this 
complex field. 

Acknowledgements. We thank the referee for his careful reading of the paper 
and his suggestions. E. H. acknowledges the support of the National Science 
Foundation (US) for the support of his program in astrochemistry, the support 
of the NSF Center for the Chemistry of the Universe, and the support of NASA 
(NAI) for his work in the evolution of pre-planetary matter. V. W. and S . G. thank 
the French program PCMI for partial support of this work. 



References 

Bergin, E. A., Hartmann, L. W., Raymond, J. C, & Ballesteros-Paredes, J. 2004, 
ApJ, 612, 921 

CardeUi, J. A., Mathis, J. S., Ebbets, D. C, & Savage, B. D. 1993, ApJ Lett., 
402, L17 

Caselli, P., Walmsley, C. M., Terzieva, R., & Herbst, E. 1998, ApJ, 499, 234 

Dickens, J. E., Irvine, W. M., Snell, R. L., et al. 2000, ApJ, 542, 870 

Flower, D. R. & Pineau des Forets, G. 2003, MNRAS, 343, 390 

Graedel, T. E., Langer, W. D., & Frerking, M. A. 1982, ApJS, 48, 321 

Indriolo, N., Fields, B. D., & McCall, B. J. 2009, ApJ, 694, 257 

Jenkins, E. B. 2004, in Origin and Evolution of the Elements, 336 — H 

Krco, M., Goldsmith, P R, & Brown, R. L. 2009, in Bulletin of the American 

Astronomical Society, Vol. 41, Bulletin of the American Astronomical 

Society, 457-+ 
Le Bourlot, J. 1991, A&A, 242, 235 
Le Petit, F, Rouelf, E., & Herbst, E. 2004, A&A, 417, 993 
Lee, H.-H., Herbst, E., Pineau des Forets, G., Rouelf, E., & Le Bourlot, J. 1996, 

A&A, 311, 690 

Lee, H.-H., Rouelf, E., Pineau des Forets, G., et al. 1998, A&A, 334, 1047 
Lepp, S. & Dalgamo, A. 1996, A&A, 306, L21 

McCall, B. J., Huneycutt, A. J., Saykally, R. J., et al. 2003, Nat, 422, 500 

Meyer, D. M., Jura, M., & Cardelli, J. A. 1998, ApJ, 493, 222 

Ohishi, M., Irvine, W. M., & Kaifu, N. 1992, in lAU Symp. 150: Astrochemistry 

of Cosmic Phenomena, Vol. 150, 171 
Osterbrock, D. E., Tran, H. D., & Veilleux, S. 1992, ApJ, 389, 305 
Padovani, M., GalU, D., & Glassgold, A. E. 2009, A&A, 501, 619 
Ruffle, D. P, Hartquist, T. W., CaseUi, P, & Williams, D. A. 1999, MNRAS, 

306, 691 

Savage, B. D. & Sembach, K. R. 1996, ARA&A, 34, 279 

Shalabiea, O. M. & Greenberg, J. M. 1995, A&A, 296, 779 

Spitzer, L. J. & Tomasko, M. G. 1968, ApJ, 152, 971 

Vasyunin, A. 1., Semenov, D., Henning, T., et al. 2008, ApJ, 672, 629 

Vasyunin, A. I., Sobolev, A. M., Wiebe, D. S., & Semenov, D. A. 2004, 

Astronomy Letters, 30, 566 
Wakelam, V. & Herbst, E. 2008, ApJ, 680, 371 
Wakelam, V, Herbst, E., & Selsis, F 2006, A&A, 451, 551 
Wakelam, V, Loison, J.-C, Herbst, E., et al. 2009, A&A, 495, 513 
Wakelam, V, Selsis, R, Herbst, E., & Caselli, P 2005, A&A, 444, 883 
Walsh, C, Harada, N., Herbst, E., & Millar, T. J. 2009, ApJ, 700, 752 
Webber, W. R. 1998, ApJ, 506, 329 

Wootten, A., Loren, R. B., & Snell, R. L. 1982, ApJ, 255, 160 



